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Abstract — Long-range  laser  communications  links  offer 
many  potential  advantages  to  satellite  communications 
systems  including  high  bandwidth,  security,  and  a  small 
physical  footprint.  However,  the  establishment  of  a  laser 
communications  link  can  be  very  difficult  and  time 
consuming  due  to  the  extremely  small  beam  divergences 
and  stringent  platform  stability  requirements  inherent  to 
lasercom  systems.  Most  of  the  difficulty  lies  in  the  first  step 
of  actually  painting  the  partner  terminal  with  a  laser  spot  or 
detecting  said  spot  due  to  large  uncertainties,  terminal 
vibrations  and  jitter,  bright  optical  backgrounds  and  clutter. 
A  detailed  physical  and  statistical  model  of  this  first 
detection  process  has  been  developed  and  verified  against 
computer  simulations  using  representative  background 
images  to  ensure  the  success  of  this  critical  first  step  in  the 
acquisition  process,  enabling  future  designers  to  have  some 
confidence  that  their  terminals  will  be  able  to  acquire12. 
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1.  Introduction 

Future  space-based  communications  systems  are 
investigating  the  use  of  laser-based  systems  due  to  the 
revolutionary  increases  in  communications  bandwidth, 
security,  and  SWaP  enabled  by  the  use  of  laser  light  as  the 
information  carrier.  Several  current  and  future  space-based 
systems,  such  as  the  ESA’s  Artemis,  the  JAXA’s  OICETs 
and  the  DOD’s  TSAT  [1-6],  are  investigating  or  have 
demonstrated  lasercom  links  between  a  satellite  and  another 
satellite  or  a  high-flying  aircraft. 


Establishing  a  lasercom  link  between  two  terminals  requires 
an  acquisition  process  to  enable  the  terminals  to  find  and 
lock  onto  each  other  with  very  small  divergence  beams  (1  to 
100  prods,  depending  on  range  and  aperture  sizes).  The 
initial  step  of  this  acquisition  process  involves  each  terminal 
scanning  an  illumination  beam  towards  its  partner  terminal 
while  trying  to  detect  and  discern  the  incoming  optical 
signal  in  the  presence  of  platform  vibrations  and  cluttered 
optical  backgrounds.  This  initial  step  of  scan  and  detect  can 
be  the  most  difficult  and  time  consuming  step  in  the  process 
of  setting  up  a  lasercom  link,  especially  when  the  system  is 
plagued  by  false  alarms  or  missed  detections. 

Early  free-space  laser  communications  systems  typically 
used  a  quad-cell  detector  for  both  acquisition  and  tracking 
[7,8],  possibly  utilizing  a  modulated  signal  to  cut  down  on 
the  large  background  noise  inherent  in  such  a  system  with 
the  field  of  view  required  to  encompass  the  uncertainty  in 
the  location  of  the  far  terminal.  As  sensor  technology  has 
evolved,  designers  have  been  drawn  towards  utilizing  CCD- 
type  focal  plane  arrays  (FPAs)  because  of  their  excellent 
responsivity,  low  noise  characteristics,  and  the  smaller 
background  signals  created  by  segmenting  the  field  of  view 
into  hundreds  to  millions  of  pixels. 

The  remainder  of  the  paper  is  organized  as  follows:  First, 
we  present  a  model  of  a  FPA  that  includes  the  relevant  noise 
sources.  Next,  a  simplified  model  to  represent  spacecraft 
jitter  is  presented.  Then  the  statistical  basis  underlying  the 
first  detection  process  is  developed,  and  after  that  the  image 
processing  steps  to  achieve  first  detection  are  developed. 
Finally,  the  example  background  image  and  simulation 
results  are  presented,  followed  by  the  conclusion. 

2.  FPA  Modeling 

The  focal  plane  array  devices  likely  to  be  used  in  laser 
communications  systems  are  Si,  InGaAs  or  other  types  of 
arrays  that  act  as  power  integrators,  converting  the  incoming 
photons  into  electrons  and  then  counting  the  number  of 
electrons  detected.  The  principles  of  operation  and  typical 
noise  sources  for  these  devices  are  well  understood  due  to 
their  extensive  use  in  astronomical  application  (see  for 
example  [9]).  The  number  of  electrons  produced  for  a  given 
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energy  input  onto  the  FPA  can  be  determined  by  the 
following  equation: 


Npe=sin-QE.^~  (1) 

where: 

Npe  =  Number  of  photoelectrons  produced 
Sin  =  Energy  input  in  1  integration  period 
QE  =  FPA  quantum  efficiency 
X  =  Wavelength 

h  =  Planck’s  constant  (6.626E-34  Js) 
c  =  Speed  of  light  (2.99792458E8  m/s) 

The  electrons  are  then  read  out  by  being  digitized  into  N 
bits,  perhaps  after  an  amplification  stage.  The  read  out  from 
the  FPA  then  becomes: 


fullwell/ 2n 

Where: 

R  =  Readout  from  the  FPA 

Fullwell  =  saturation  level  of  a  FPA  pixel,  in 

electrons 

n  =  Number  of  bits  in  digitizer 


changes  in  the  background  image  itself  It  can  be  modeled 
as  some  percentage  of  the  background  signal,  where  cf  is  the 
clutter  factor: 

Nbc  =cf*  background  (4) 


Digitization  noise  arises  from  binning  the  electron  counts  in 
the  analog  to  digital  converter,  and  is  given  by: 


Nad  = 


Full  Well 


(5) 


Pixel  non-uniformity  noise  arises  from  the  fixed  pattern 
pixel-to-pixel  variations  in  responsivity.  This  can  be  quite 
large  initially;  however,  careful  calibration  with  a  flat  field 
can  reduce  it  to  a  small  value.  Unfortunately,  the  spatial 
nonuniformity  is  not  constant  and  can  drift  slowly  in  time 
(see  [16]).  The  noise  term  can  be  given  in  terms  of  a 
nonuniformity  factor  nu: 

Nnu  =  (nu  *  ( signal  +  background  +  dark))  (6) 

The  total  noise  per  pixel  is  then  given  by: 

^ = ENNNNNzNkNnl  (7) 


Several  sources  of  noise  conspire  to  make  life  more 
difficult.  Typical  noise  sources  are  described  in  detail 
below  and  include  shot  noise,  read  out  noise,  dark  current 
noise,  background  clutter  noise,  digitization  noise,  and  pixel 
non-uniformity  noise. 

Shot  noise  is  a  Poisson  process  that  arises  from  the  quantum 
mechanical  nature  of  photon  arrival  counting,  and  can  be 
represented  as: 

Nshot  =  yl signal  +  background  +  dark  (3) 


Where: 

Signal  =  #  of  signal  electons 

Background  =  #  of  background  image  electrons 

Dark  =  #  of  dark  current  electrons 


Computing  the  SNR  for  the  centroiding  window  requires 
summing  up  the  signal  and  noise  terms  in  each  pixel,  and  is 
given  by: 


SNR  = 


gA 


(8) 


An  estimate  of  the  Noise  Equivalent  Angle  (NEA)  is 
proportional  to  the  inverse  of  the  SNR: 


NEA  = 


SNR 


ifov 


(9) 


where  rj  is  a  dimensionless  conversion  factor.  Ortiz,  Lee  and 
Hemmati  [10]  give  r|  to  be: 


Read-out  noise  is  associated  with  the  process  of  reading  out 
the  electron  well.  It  can  be  rate  dependent,  but  is  considered 
to  be  independent  of  rate  in  this  simulation.  The  value  is 
device  dependent,  but  values  of  tens  to  a  few  thousand 
electrons  are  typical. 

Dark  current  noise  is  caused  by  the  thermal  energy  of  the 
FPA.  It  becomes  important  as  the  wavelength  is  increased 
and  integration  times  grow  longer.  For  the  wavelengths  and 
short  integration  times  inherent  to  lasercom  systems,  it  can 
be  assumed  to  be  negligible. 

Background  clutter  noise  is  the  residual  image  that  is  left 
over  after  the  background  has  been  subtracted  from  the 
image,  and  can  be  considered  as  arising  from  three  sources: 
errors  in  registering  the  background  image;  motion  of  the 
optical  line  of  sight  with  respect  to  the  background;  and 


1  - 


vAo+i) 


(10) 


where  N0  is  the  truncated  half  width  of  the  window.;  but  we 
believe  this  may  be  a  misprint  and  propose  an  estimation  of 


rj  as: 


71  = 


3 

*>„+i) 


(ii) 


where  Np  is  the  window  size;  However,  results  from  our 
simulations  show  this  to  be  a  very  good  estimator  only  for 
some  combinations  of  window  size  and  SNR;  in  other  cases 
it  can  be  off  by  up  to  a  factor  of  3.  This  form  of  the  NEA 
equation  illustrates  the  trade  that  must  be  made  to  determine 
the  window  size:  as  Np  goes  up,  r|  decreases,  leading  to  a 
smaller  NEA,  but  the  SNR  also  drops  with  increasing  Np, 
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(13) 


which  leads  to  a  larger  NEA.  The  optimal  window  size  can 
be  determined  based  on  an  analysis  of  the  expected  signal 
and  noise  terms 


oo 


Pd  =  \P,nMYlx 

th 


3.  Jitter  Modeling 


Optical  communications  payloads  are  typically  secondary 
payloads  on  larger  spacecraft  and  are  subject  to  the  jitter 
vibration  environment  caused  by  other  onboard  devices  such 
as  momentum  wheels,  solar  arrays,  thrusters,  etc.  In  the 
design  process,  the  actual  or  expected  vibration  spectrum  for 
the  platform  the  lasercom  payload  is  destined  for  should  be 
used.  Unfortunately,  only  a  few  data  points  containing 
satellite  rotational  vibration  data  are  available  [11,12].  For 
the  purposes  of  this  paper,  a  typical  vibration  environment 
from  [7]  was  assumed,  with  the  PSD  shown  in  the  top  half 
of  Figure  1: 


Figure  1:  Jitter  PSD  (1-axis)  and  Time  History 


A  time  history  was  created  from  this  PSD  by  applying  an 
appropriately  weighted  inverse  Fourier  transform  to  the 
PSD,  with  a  random  phase  assigned  to  each  frequency,  as 
shown  in  the  bottom  half  of  Figure  1.  Both  X  and  Y  time 
histories  were  created,  and  were  uncorrelated  when  used  in 
the  simulations  later  in  the  paper. 

4.  Detection  Statistics 


Simple  detection  theory  as  described  in  [13,14]  is  used  to 
build  up  the  detection  statistics  for  the  entire  FPA.  The  basis 
for  this  is  shown  in  Figure  2  and  depends  on  the  PDFs  for 
the  noise  and  the  signal  (Pn  and  Pn+S),  the  probability  of  false 
alarm  (pfa),  the  probability  of  detection  (pd)  and  the 
threshold  level  (th).  For  a  given  noise  PDF,  setting  the 
threshold  level  determines  the  probability  of  false  alarm  by: 


oo 

Pfa=\Pn{x)dx 


(12) 


Similarly,  the  probability  of  detection  is  determined  by: 


Equations  (12,13)  are  for  a  single  detector  at  a  single  point 
in  time  (designated  as  Pfa.p,  probability  of  false  alarm  per 
pixel).  When  using  a  FPA  as  a  detector,  the  designer  is 
typically  given  a  higher  level  system  specification  such  as 
“a  1%  chance  of  a  false  alarm  per  minute”  (Pfa.s,  probability 
of  false  alarm,  system),  which  requires  computing  the  false 
alarm  rate  over  all  pixels  and  a  number  of  frames.  The 
required  probability  of  false  alarm  per  frame  (Pfa.f)  can  be 
derived  from  the  system  probability  by: 


fa-s 

FRT 


(14) 


where  FR  is  the  frame  rate  and  T  the  time  interval  for  the 
specification  (i.e.,  one  minute).  Computing  the  required 
pixel  false  alarm  rate  from  the  frame  false  alarm  rate  is 
accomplished  by: 

Pfa-f  Pfa-PYP  (15) 


Where  Np  is  the  number  of  pixels  per  frame,  N2.  This 
assumes  a  false  alarm  is  declared  for  a  given  frame  even  if  a 
true  detection  occurs  in  some  other  pixel  in  the  frame. 
Combining  Equations  (14,15)  and  solving  for  Pfa.p  gives: 


pp-P=  i- 


i— 


fa-s 

FR-T 


1/np 


(16) 


As  the  specification  for  the  required  false  alarm  rate 
becomes  more  restrictive,  the  required  threshold  setting  will 
move  to  the  right  in  Figure  2.  The  requirement  for  the 
probability  of  detection  (Pd),  coupled  with  the  level  the 
threshold  is  set  at,  determines  how  far  to  the  right  the  PDF 
Pn+S  must  shift,  which  is  realized  by  setting  a  requirement  on 
the  received  signal  to  noise  ratio  (SNR). 


3 


Most  systems  will  have  a  limit  to  the  achievable  SNR, 
especially  in  the  beginning  stage  of  acquisition,  which  by 
the  above  analysis  sets  a  limit  on  the  false  alarm  rate  that 
can  be  realized  for  any  given  sensor.  Fortunately,  there  are 
other  design  options  involving  multiple  pixels,  frames,  and 
signal  patterns  that  can  be  exploited  to  reduce  the  required 
SNR  while  still  meeting  requirements  for  Pfa  and  Pd.  For 
example,  if  the  signal  is  required  to  be  present  for  nf  frames 
(consecutive  or  in  some  predetermined  sequence),  Equation 
(15)  becomes: 

Pfa-f=^-PnflJP  (17) 

Alternatively,  you  can  utilize  the  optical  design  of  your 
sensor  to  create  a  point  spread  function  that  will  encompass 
multiple  pixels.  Of  course,  in  this  case  the  trade  is  a  lower 
SNR  in  each  pixel.  However,  you  are  likely  to  require  the 
PSF  to  span  multiple  pixels  in  any  case  for  centroiding 
accuracy  reasons.  For  nc  contiguous  pixels,  Equation  (15) 
becomes: 

Pfa_f=\-(\-p;;_pf"‘  (18) 

In  this  case,  the  number  of  pixels  Np  is  replaced  by  the 
number  of  nc  contiguous  pixel  groupings  Npnc.  For  nc  =  2: 

Nnn  =  4N2  -6A  +  2  (19) 

Pnc 

To  illustrate  the  advantage  schemes  such  as  these  provide, 
the  following  situation  is  analyzed: 

N  =512  (number  of  pixels  in  1  direction) 

FR  =  60  Hz  (frame  rate  of  camera  in  acquisition 
mode) 

T  =  1 80  seconds  (time  over  which  acquisition  is 
attempted) 

Pfa-s  =  1%  (required  system-level  false  alarm 
probability) 

Using  these  parameters,  the  values  in 
Table  1  were  generated  for  nc  =  1,2,3  and  nf  =  1,2,3.  It  is 
apparent  that  vast  reductions  in  the  required  false  alarm  rate 
(and  hence,  SNR)  can  be  achieved  by  these  techniques. 


Per  pixel  false  alarm  requirements 

#  of 

1 

#  of  frames 

2 

3 

pixels 

1 

3.50E-12 

1 .90E-06 

1 .50E-04 

2 

9.40E-08 

9.70E-04 

9.80E-03 

3 

9.60E-05 

9.80E-03 

4.60E-02 

Table  1:  Derived  False  Alarm  Rates,  Per  Pixel 


5.  Image  Synthesis  and  Processing 


Several  image  processing  steps  are  required  to  go  from  the 
input  image  to  a  time  sequence  of  noisy,  jittering  images, 
and  then  to  perform  the  background  suppression  and 
acquisition  detection  functions.  For  the  purpose  of  the 
acquisition  simulation  developed  here,  the  input  image  is 
assumed  to  be  a  ‘truth’  image,  devoid  of  any  noise  effects. 
Jitter  and  noise  are  then  added  to  the  input  image  to  produce 
a  sequence  of  input  images,  and  a  laser  spot  with  the  desired 
characteristics  is  synthesized  and  added  to  the  sequence  at  a 
random  time  stamp.  Finally,  the  image  sequence  is  analyzed 
to  detect  an  acquisition  by  using  background  suppression 
and  the  detection  algorithms  discussed  in  the  previous 
section. 

Jitter  is  added  to  the  image  sequence  by  taking  the  jitter  time 
histories  as  developed  in  Section  3.  The  jitter  files  are  first 
converted  from  microradians  to  pixels  by  dividing  by  the 
pixel’s  instantaneous  field  of  view  (iFOV).  A  several  pixel 
wide  buffer  zone  is  removed  from  the  edges  of  the  input 
image  to  allow  for  translation  and  drift.  If  the  original  image 
is  denoted  as  A,  and  a  translation  of  (X,Y)  whole  pixels  is 
denoted  as  AXjy,  while  a  translation  of  partial  pixels  is 
denoted  as  Ay.  To  produce  the  image  that  is  shifted  to 
Ax+iy+j,  the  image  is  first  shifted  by  (x,y)  pixels  and  then 
appropriately  interpolated  with  the  image  shifted  to 
Aciei(x+i);Ciei(y+j) . 

The  more  difficult  task  is  to  add  in  noise  to  each  frame. 
Each  shifted  frame  is  used  as  a  basis  for  adding  in  per  pixel 
noise  based  on  the  noise  processes  described  in  Section  2. 
The  standard  deviation  of  the  noise  an  from  Equation  (7)  is 
used  to  randomly  generate  a  new  pixel  value  from  Equation 
(1)  as  follows: 


Npe{k)=v(k) 


(20) 


Where: 


E\y{k)\  =  sin  -Qe'T-’  var[v(&)]  =  cr2  (2 1 ) 

nc 


Once  a  sequence  of  noisy,  jittered  images  has  been 
generated,  a  simulated  laser  spot  is  generated  using  an  Airy 
profile  as  given  by: 


f  r 

2  Jx 


m=ic 


A 


nOD 


eff 


v  v 


(22) 


z  z 


With: 


I=- 


PttD 


eff 


4A1 


(23) 


Where  P  is  the  input  power  and  Deff  determines  the  radius  of 
the  Airy  spot.  The  spot  is  randomly  located  within  the  array 
at  the  point  (Cx,Cy)  and  an  integration  is  performed  over  the 
boundaries  of  each  pixel  as  follows: 
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x  max,-  y  max j  ,  - - - - — x 

S,.J=  \  \  lU(x-Cxy+(y-Cj)dydx  (24) 

xmirij  y  min  y  v  ' 

The  sequence  of  frames  is  analyzed  by  first  removing  the 
strong  background  signal  from  the  images.  Several 
techniques  were  analyzed,  depending  on  which  detection 
algorithm  was  chosen,  but  all  are  of  the  form: 

Bt=At-At~m  (25) 

Detection  is  then  simply  based  on  threshold  setting  and 
applying  the  specified  algorithm  (ie,  2  consecutive  frames  or 
2  contiguous  pixels). 

6.  Example  and  Simulation  Results 

An  example  image  was  used  as  a  basis  for  a  simulation  of 
this  process  from  the  Enhanced  Thematic  Mapper  Plus 
(ETM+)  camera  on  Landsat  7  [15].  This  camera 

encompasses  seven  spectral  bands,  ranging  from  .45 
micrometers  to  2.35  micrometers;  an  image  from  band  2 
(.52  to  .6  p)  was  chosen.  Also,  for  this  satellite,  the 
Instantaneous  Field  of  View  (IFOV)  equals  42.5 
microradians  equating  to  a  nominal  ground  resolution  of  30 
meters.  Landsat  7  satellite  has  a  40.64  cm  aperture  size  with 
a  243.8  cm  focal  length.  The  example  image  is  shown  in 
Figure  3: 


Figure  3:  Example  Reference  Image 


This  image  has  been  optimized  for  imaging  purposes,  as  the 
histogram  shown  in  Figure  4  clearly  shows  (the  background 
brightness  values  cover  the  full  range  of  the  camera).  A 
sensor  with  parameters  optimized  to  be  a  lasercom  detector 
might  produce  images  that  look  much  darker;  it  is  unlikely 


that  you  would  want  to  increase  the  brightness  as  it  becomes 
impossible  to  detect  an  additional  signal  in  a  saturated  pixel. 
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Figure  4:  Image  Histogram 

A  sequence  of  frames  was  developed  by  applying  the 
vibration  data  described  in  Section  3  in  the  manner 
explained  in  Section  5.  An  example  of  the  shifting  that 
occurs  is  shown  in  Figure  5,  where  the  motion  has  been 
greatly  magnified  for  effect.  These  frame  to  frame  shifts  are 
a  major  source  of  noise  (background  clutter)  when  trying  to 
detect  a  weak  signal  via  background  subtraction  techniques. 


Figure  5:  Image  Motion  Sequence 
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Noise  was  added  on  a  per-pixel,  per  frame  basis  as 
described  in  Section  5,  with  the  following  input  parameters: 
Read  noise:  39  electrons,  1-g 
Pixel  non-uniformity:  0.02% 

A/D  conversion:  8  bits 

An  example  with  5x  the  actual  noise  added  is  shown  in 
Figure  7,  which  is  a  blow-up  of  the  top  left  portion  of  the 
image  in  Figure  3  (the  corresponding  noise-free  image  is 
shown  in  Figure  7).  Without  the  5x  amplification  in  noise,  it 
is  very  difficult  to  visually  detect  the  degradation  in  the 
image.  It  should  be  noted  that  the  dominant  noise  source  for 
these  parameters  is  signal  shot  noise. 


Figure  6:  Image  with  Added  Noise 


Figure  7:  Image  Without  Added  Noise 


A  sequence  of  900  frames  (5  seconds  at  180  Hz)  was 
created  by  applying  the  vibration  spectra  and  random  noise 
inputs  as  described,  and  then  a  change  file  was  created  by 
subtracting  off  the  previous  image  from  each  successive 
image.  On  any  given  image,  the  histogram  of  the  remaining 


noise  in  the  change  file  looks  like  Figure  8,  while  the 
standard  deviations  in  the  sequence  of  residual  noises  are 
shown  in  Figure  9;  the  standard  deviation  ranged  from  a 
minimum  of  3.5  to  a  maximum  of  35.8,  with  the  average 
standard  deviation  being  1 1.6. 


Figure  8:  Background  Subtracted  Histogram 


35 


°0  100  300  300  <00  500  600  TOO  800  900 

Frame 

Figure  9:  Background  Subtracted  Noise  SD  per  Frame 

Once  statistics  on  the  residual  noise  are  known,  it  then 
becomes  possible  to  determine  the  required  signal  levels  to 
achieve  a  given  false  alarm  rate  and  probability  of  detection. 
Let  us  suppose  a  1%  false  alarm  rate  over  our  5  second 
window  is  required.  Using  Equation  (15)  gives  a  per- frame 
false  alarm  rate  of  1.11E-5.  With  a  simplistic  1 -pixel,  1- 
frame  detection  algorithm,  the  required  per-pixel  false  alarm 
rate  is  3.9E-11  from  Equation  (16),  which  requires  a 
threshold  set  to  a  level  greater  than  6.5a.  This  equates  to  a 
threshold  at  75  if  the  mean  standard  deviation  is  used  and 
233  if  the  maximum  value  is  used.  The  maximum  pixel 
value  per  frame  is  shown  in  Figure  10;  note  that  a  threshold 
of  75  will  produce  a  plethora  of  false  alarms,  while  a  level 
of  233  produces  none  where  we  would  expect  to  see  one 
only  1%  of  the  time.  Using  neither  the  mean  nor  the 
maximum  standard  deviation  has  achieved  the  desired 
results,  so  a  closer  look  at  what  is  going  on  in  Figure  9  and 
Figure  10  is  warranted. 
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Figure  10:  Maximum  Change  in  Pixel  Values  per  Frame 

The  (cropped)  image  we  are  examining  contains  531 2 
pixels;  therefore  it  should  not  be  surprising  that  events  with 
probabilities  as  low  as  3.5E-6  occur,  which  corresponds  to 
almost  4.5  a.  However,  Figure  11  shows  the  number  of 
standard  deviations  the  maximum  value  is  away  from  the 
mean  in  each  frame;  events  are  occurring  at  over  12  a  away 
from  the  mean,  which  corresponds  to  a  probability  of  <  1E- 
32.  Obviously  what  is  going  on  is  that  the  distribution 
shown  in  Figure  8  is  not  strictly  normal  in  the  tails. 
Examination  of  where  the  extremes  occur  in  the  image 
shows  that  sharp  boundaries  between  high  and  low  valued 
pixels  are  the  main  sources  of  frame-to-frame  variation.  In 
any  case,  for  this  particular  algorithm,  the  false-alarm 
statistics  built  up  in  Section  4  can  not  be  used  in  this  manner 
to  produce  a  valid  false  alarm  rate  based  threshold.  Instead, 
one  would  have  to  take  the  data  shown  in  Figure  10  and  fit 
it  to  an  extreme  value  distribution  to  determine  the  correct 
threshold. 


Frame 


The  next  logical  step  to  take  is  to  require  the  spot  to  remain 
for  2  frames.  From  Equation  (18),  this  gives  a  per-pixel 
false  alarm  rate  of  6.28E-6,  which  corresponds  to  a 
threshold  set  to  4.4a.  Changing  the  detection  criteria  also 
changes  the  residual  noise  statistics  (the  a)  from  the  single 
frame  detection  case,  in  a  manner  that  is  dependent  on  how 
well  the  residual  noise  is  correlated  from  frame  to  frame  and 
how  the  background  subtraction  is  performed.  For  example, 
if  the  scheme  is  =  An  -  An_1,  B2n  =  An  -  An"2,  then  the 
noise  in  B2n  will  tend  to  be  larger,  as  more  jitter  has 
occurred.  If  a  new  frame  Cn  is  constructed  by  finding  the 
minimum  value  for  each  pixel  in  Bin,  B2n,  the  standard 
deviation  over  the  run  time  becomes  as  shown  in  Figure  12, 
which  is  only  moderately  reduced  from  Figure  9;  it  has  a 
min  of  3.2,  a  max  of  33.3,  and  an  average  of  12.11;  the 
maximum  values  in  Cn  are  likewise  only  slightly  reduced 
from  those  shown  in  Figure  12.  This  indicates  that  this 
particular  scheme  produces  noise  in  successive  frames  that 
is  highly  correlated,  and  the  assumptions  used  to  come  up 
with  the  4.4a  threshold  level  are  suspect  (Setting  the  level  to 
4.4a  when  a  is  the  max  value  of  33.3  gives  a  threshold  at 
147;  this  produces  100  false  alarms  during  this  period!). 


Figure  12:  Two  Frame  Noise  SD  per  Frame 


An  alternate  scheme  to  achieve  the  desired  two  frames  false 
alarm  rates  is  to  let  B,n  =  A"  -  A"'1,  B2n  =  A"1  -  A"'2.  In  this 
case,  the  noise  becomes  much  less  correlated  and  the 
standard  deviations  are  reduced  to  the  values  shown  in 
Figure  13,  with  a  min  of  2.7,  a  max  of  22.8,  and  an  average 
value  of  9.8.  The  maximum  pixel  value  that  is  present  in 
both  frames  is  shown  in  Figure  14;  it  is  clear  that  a  much 
lower  threshold  can  be  tolerated  in  this  case. 


Figure  11:  Maximum  value/Standard  Deviation 
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Figure  13:  Alternate  Two  Frame  SD 


Figure  14:  Alternate  Two  Frame  Max  Change  per 
Frame 


Setting  the  threshold  at  4.4a,  and  using  the  maximum  a  of 
22.8,  gives  a  level  of  97;  the  maximum  value  of  the  noise 
does  not  exceed  this  level  in  any  of  the  frames. 

Finally,  combing  the  multi-frame  criteria  shown  above  with 
the  pair  of  pixels  criteria  as  described  in  Equations  18-19 
reduces  the  observed  maximum  values  even  more.  The 
results  for  this  simulation  are  shown  in 


Figure  16.  Note  that  the  maximum  values  have  been 
reduced  by  about  10%  from  the  previous  case  considered. 
The  per-pixel  false  alarm  rate  has  been  dropped  all  the  way 
to  1.8E-3,  which  corresponds  to  a  threshold  set  to  2.9a.  The 
standard  deviation  for  this  process,  shown  in 


Figure  16,  has  a  mean  value  of  10.8  and  a  maximum  of 
27.6.  As  the  threshold  has  been  set  at  2.9a,  the  non-normal 
behavior  in  the  tails  of  the  distribution  is  not  apparent,  and 
we  are  much  closer  to  what  is  expected:  At  a  threshold  of 
2.9a  where  the  maximum  sigma  is  used,  we  get  5  instances 
of  threshold  exceedance,  which  is  more  than  we  would 
expect. 
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7.  Conclusions 


Figure  15:  Two  Frame,  Two  Pixel  Max  Values 


Frame 

Figure  16:  Two  Frame,  Two  Pixel  SD 

Table  2  summarizes  the  results  from  this  section.  The  first 
column  lists  each  algorithm  that  was  tested,  along  with  the 
required  per-pixel  false  alarm  rates.  Type  A  corresponds  to 
the  single-frame  algorithm,  Type  B  the  first  two-frame 
algorithm,  Type  C  the  second  two-frame  algorithm,  and 
Type  D  the  two-frame,  two-pixel  algorithm.  The  second 
column  indicates  the  estimated  threshold  level  to  achieve 
the  desired  false  alarm  rate  based  on  the  mean  standard 
deviations,  while  the  third  column  is  based  on  the  maximum 
recorded  standard  deviation.  The  final  column  represents  the 
maximum  value  recorded  in  each  simulation. 


Threshold  Setting  Results 


Threshold  levels  (1%  FAR) 


Type,  PPFAR 

Est  (mean) 

Est  (max) 

Actual 

A  /  3.9E-11 

75 

233 

231 

B  /  6.28E-6 

53 

147 

208 

C  /  6.28E-6 

43 

100 

101 

D/1.8E-3 

31 

80 

88 

An  acquisition  detection  simulation  was  developed  for 
detecting  faint  point  source-like  signals  against  bright  and 
highly  cluttered  background  environments.  The  simulation 
incorporated  electronic  and  other  inherent  noise  sources, 
and  accounted  for  vibration-induced  image  jitter.  Detection 
statistics  were  developed  for  several  different  detection 
algorithms,  and  an  example  using  a  representative  Landsat 
image  was  presented  to  demonstrate  the  effectiveness  of  the 
simulation  and  the  derived  detection  statistics. 

The  simulation  found  that  observed  maximum  threshold 
values  fell  as  expected  with  decreasing  per-pixel  false  alarm 
rate  requirements  as  more  restrictive  detection  algorithms 
were  used.  However,  it  was  found  that  the  equations 
developed  in  Section  4  did  not  provide  an  absolute 
quantitative  match  to  determining  the  expected  required 
threshold  levels.  This  was  found  to  be  primarily  due  to  non¬ 
normal  behavior  in  the  tails  of  the  distributions,  caused  by 
correlations  in  bright  and  dark  pixels. 

Future  work  will  involve  modifying  the  simulation  for 
tracking  applications  as  well  as  detecting  multiple  sources 
when  very  bright  sources  are  present.  Further  development 
of  the  detection  statistics  to  more  accurately  match  the 
observed  results  will  also  be  performed. 
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